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ABSTRACT 

Generically, massive gravity gives a non-unique gravitational field around a star. For a special 
family of massive gravity theories, we show that the stellar gravitational field is unique and 
observationally acceptable, that is close to Einsteinian. The black hole solutions in this family of 
theories are also studied and shown to be peculiar. Black holes have a near-horizon throat and 
the curvature diverging at the horizon. 

We show that there exists a sub-family of these massive gravity theories with non-singular at 
horizon black holes. 



1. Introduction 

The history of massive gravity, as it relates to 
the present work, can be summarized as follows. 

Fierz and Pauli (1939) introduce linearized 
massive gravity. 

Van Dam and Veltman (1970) and Zakharov 
(1970) show that the Fierz-Pauli linearized mas- 
sive gravity is strongly ruled out by the solar sys- 
tem (and many other) observations - the so-called 
vDVZ discontinuity. 

Vainshtein (1972) argues that, due to non-linear 
effects, massive gravity with small graviton mass 
is in fact observationally acceptable - it gives 
close to Einsteinian gravitational field near mas- 
sive enough bodies. 

The Vainshtein (1972) argument is not a real 
proof. Vainshtein just points out that a near- 
Einsteinian solution of a spherically symmetric 
static problem exists in a range of radii. But the 
existence of an acceptable solution in a finite range 
of radii does not, generically, guarantee that the 
same solution has the right asymptotic behavior 
at large and small radii. And indeed, Babichev et. 
al. (2010) show that the spherically symmetric 
solution in massive gravity is not even unique. 

A non-unique gravitational field of a massive 
body is strange and might be taken as a sign that 
Fierz-Paulian massive gravities (all possible non- 
linear completions of the original Fierz-Pauli ac- 



tion, call it FP) are sick. In fact ambiguous stellar 
gravity is not even the main problem with FP. FP 
is known to have ghosts (Boulware, Deser 1972) 
and unacceptable cosmologies (Gabadadze, Gruzi- 
nov 2005). 

It is not clear how general are the no-go results 
of Boulware and Deser (1972) and of Gabadadze 
and Gruzinov (2005). Maybe one can find some 
very special FP theory which is free of ghosts and 
has an acceptable cosmology. An important re- 
cent finding (De Rham et al 2010) is that there 
exists just a two-parameter sub-family of FP the- 
ories which might be ghost free - call it FP2 (Ap- 
pendix A). All other FP theories have ghosts. 

In this paper we show that 

• FP2 gives a unique spherically symmetrical 
gravitational field for a given star or black 
hole. 

• Some FP2 theories give observationally ac- 
ceptable gravitational field for stars. 

• Black holes in generic FP2 are peculiar - 
they have a throat and a singular horizon. 

• There is a one-parametric sub-family of FP2 
- call it FPl - which gives black holes 
with non-singular at horizon (and close to 
Schwarzschild for small graviton mass) ge- 
ometry, although there still exists a scalar 
quantity which is singular at horizon. 
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While the status of both FP2 and FPl is almost 
as questionable as that of FP (being FP2 is a nec- 
essary condition for ghost-freedom, but it is not 
known if it is also sufficient), we present our re- 
sults hoping that they might be qualitatively use- 
ful for any future theory of massive gravity. If the 
FP2 family does contain a good theory of massive 
gravity, then we tentatively predict that this good 
theory is from the FPl sub- family. 

2. Uniqueness of the static spherical solu- 
tions in FP2 

Representing the contribution of the FP action 
terms by an additional stress-energy tensor in the 
Einstein equation, we write the spherically sym- 
metric solution in the standard form 

ds^ = e^dt^ - e^dr^ -r^dVi^. (1) 

One has to introduce an additional unknown 
(Vainshtein's ^(r)) to describe the FP stress- 
energy. In FP2, however, the resulting system of 
equations is of lower order than in the generic FP. 
As shown in Appendix B, the vacuum equations 
can be written as 

v'^h{r,v,\), (2) 

A' = /2(r,i^,A), (3) 

where the functions /i , /2 are given (implicitly) in 
Appendix B. 

Equations p|3p require two integration con- 
stants, in contrast to GR where the A-equation 
decouples and only one integration constant is 
needed. However, the condition of asymptotic flat- 
ness (z^ — >■ 0, A — )■ for r oo) fixes one of the 
integration constants. Then, just like in Einstein 
theory, one gets a one-parameter family of vacuum 
solutions. This remaining parameter (mass in the 
Einstein theory) characterizes the central object. 

3. Stars in FP2 

At r ^ m~^, we start with the linearized 
asymptotically flat solution (given in Appendix 
B) and numerically integrate the system ^ ^ 
inwards. We modify the equations ^ ^ appro- 
priately after we cross the surface of the star in 
order to include the stress-energy of matter (Ap- 
pendix C). For a given stellar radius and mass (as 



determined by the large-r asymptotic), the correct 
value of the stellar energy density e is obtained by 
numerical shooting. We use a simple matter equa- 
tion of state e = const. 

If the gravity parameters are chosen appropri- 
ately 0, and so long as the graviton mass is small 
enough, the star solution agrees with the Einstein 
theory at small radii, and it agrees with the lin- 
earized massive gravity result at large radii (Fig. 
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-Xjv 

— z 

































io"3 10"^ 10"^ 10° 



Fig. 1. — Semi- logarithmic plot of —\/u and redshift 
2 = e~"l'^ — 1 as a function of r. —\/u linearly grows with 
r at large values of r in agreement with the linearized mas- 
sive gravity solution IIB17I I . Einstein theory is recovered at 
small radii since —\/v approaches the Schwarzchild metric 
value (i.e. X/v = —1) outside the star, and A — >■ at the 
center. The surface of the star is located at Ts = 0.01 and 
its mass is determined implicitly hy v = — 0.005e~"""/r at 
large r. The gravity parameters are m = 1, C2 = 0, C3 = 2. 

4. Black Holes in FP2 

For black holes, equations ^ [3]) cannot be in- 
tegrated as written, because the solution features 
a coordinate singularity w = e~^l'^ = at a finite 
redshift z = e-'^/^ _ i_ 

This means that the black hole develops a 
throat - as one moves inward to larger z, the cir- 
cumference of the sphere z = const first decreases 



As seen from the Bianchi identity 1 IBI6I I. to recover the 
linearized Schwarzschild solution, one needs to choose the 
gravity parameter C3 > 0, in agreement with the decou- 
pling limit analysis of Chkareuli and Pirtskhalava (2011). 
A careful study of the parameter space (c2,C3) was not 
performed. In what follows, C3 > and |c2| is not too big. 
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but then increases again. The solution then can- 
not be calculated with the coordinatization ([T]). 
We changed the coordinatization to 



(4) 



where p is the new independent variable, while 
u = v[p) and r = r[p) are the new unknowns. We 
re-write the system ([21 E]) as (prime now denotes 
the p-derivative) 

r' ~ w (5) 



f2{r,w,z). 



(6) 
(7) 



These equations can be integrated all the way to 
the horizon z = oo. The results are shown in Fig. 2. 
As seen, one does recover the Schwarzschild metric 
at small radii, except in the near-horizon region. 
Near the horizon, the black hole develops a throat, 
and the curvature diverges at the horizon. 

These features can be understood, in the small 
graviton mass limit, through an approximate near- 
horizon solution (Appendix D). Measuring p from 
the point of minimum circumference 27rrmin, we 
get the (circumference-defined) radius r and the 
Ricci curvature invariant W = R^^R^^'- 



W ~ Cim'^z 



n 

C2m? /{p/r^n 



2J2 



(8) 
J (9) 
(10) 



Here Ci are positive dimensionless numbers which 
depend on the dimensionless parameters of the 
theory C2 and C3. As p decreases from positive to 
negative values, the circumference reaches a min- 
imum value outside the horizon. Then one ap- 
proaches the horizon, z = 00, and W diverges. 

As shown in Appendix D, All Ci actually vanish 
for C3 = This defines FPl - a one-parameter 
sub-family of FP2. 

5. FPl 

We do not understand what makes the FPl 
(FP2 with C3 = C2, see Apendices A, B for the def- 
initions) special in any intuitive way. But we do 
understand which special property of the FPl fam- 
ily is responsible for the non-singular black hole 
horizon. 
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Fig. 2. — (1 + z)dr/dp as a function of (1 + z)~^. Here 2 
is the redshift, r is the radius as measured by the circum- 
ference of a constant-redshift sphere, p is the proper radial 
distance. For the Schwarzschild metric, (1 -|- z)dr/dp = 1, 
which is indeed seen at intermediate redshifts. The nega- 
tive values of (1 -|- z)dr/dp near the horizon {z = 00) cor- 
respond to a throat. The gravity parameters are m = 1, 
C2 = 0, C3 = 2. The mass of the black hole is given implic- 
itly hj u = — 0.06e~'"'"/r at large r. 



The FPl family is distinguished by the exis- 
tence of a unitary-gauge orthogonally-coordinatized 
de-Sitter (and Schwarzschild-de Sitter) solutions, 
Nieuwenhuizen (2011). Namely, FPl has solu- 
tions: 



ds^ = Adt^-A-^dr^-Br^{d9^+sm^ 



'% (11) 



(12) 



r 3c2 [l + C2y 

<j)^ = ~ {t^r sin cos (/>, r sin sin 0, r cos 9). 

(13) 

Here the Schwarzschild radius is arbitrary, m 
is the graviton mass, and C2 is the dimensionless 
parameter of the FPl family. 

The existence of the exact Schwarzschild-de Sit- 
ter solution is of no direct interest for the present 
study - we are only interested in the asymptoti- 
cally fiat solutions. But, as shown in Fig. 3, in FPl 
the asymptotically fiat solutions switch to the ex- 



3 



Acknowledgment s 









// \ 

* 








// 


\- 


f 




















W/4 



0.0 0.2 0.4 0.6 0.8 1.0 



(1 + 2)-' 

Fig. 3. — A2, w{l + 2) — 1, indicating closeness to 
Schwrazschild, and ^ R^yRi^'^ / i are shown vs (1 + z)~^ . 
Thin: C2 = 0.9, thick: C2 = 0.999. Other gravity parame- 
ters: m = 1, C3 = 1. The mass of the black hole is given 
implicitly hy v = — 0.04e~™'^/r at large r. 

act Schwarzschild-de Sitter solution at redshift 
. = -. (14) 

6. Conclusion 

There is a one-parameter family of Fierz- 
Paulian gravities, FPl, which might have "good" 
black holes. As shown in Fig. 3, these black holes 
are still strange. The spacetime switches to exact 
Schwarzschild-de Sitter at a finite redshift. As far 
as we can see, however, these black holes are the 
only non-singular black holes from the FP2 family 
of De Rham et al (2010). Given that FP2 is the 
only candidate ghost-free Fierz-Paulian gravity, 
these strange black holes might deserve further 
studies. 



^ The solution shown in Fig. 3 was obtained as the c2 = 
y/c^— limit of the generic FP2 solution. For exactly FPl, 
the radius corresponding to z = l/c2 becomes a bifurcation 
point, i.e. at smaller radii {z > I/C2) one can either stay on 
the Schwarzschild-de Sitter branch A2 = — l/c2, or proceed 
to even lower values of A2. The latter behavior can be 
enforced by canceling the common factor (1 -|- C2A2) from 
both sides of the Bianchi identity IIBI6I 1. This branch has 
infinite curvature at horizon. 



We thank L. Berezhiani, G. Chkareuli, S. 
Dubovsky, G. Gabadadze, S. Mukhanov, and D. 
Pirtskhalava for useful discussions. 

A. Massive gravity 

Massive gravity can be represented as Einstein 
gravity interacting with 4 non-canonical scalar 
fields (Arkani-Hamed et al 2003, Dubovsky 2004, 
Chamseddine, Mukhanov 2010). The action is the 
Einstein-Hilbert action plus the (non-linear ver- 
sion of) Fierz-Pauli action: 

S = Seh + Sfp, (Al) 
SEH^-\jd''x^R, (A2) 
SFP = m^ ^ d^x^f^U. (A3) 

Here U is an arbitrary symmetric function of the 
eigenvalues of the matrix 

- ?7i3C5^''5^0^5.0^, (A4) 

where g^i, is the metric tensor, rjAB = diag(l, — 1, 
— 1,-1) is the Minkowski matrix, and t/)'^, A = 
0, 1, 2, 3, are the 4 scalar fields. 

Since the action is generally covariant, one is 
free to choose an arbitrary coordinatization. In 
particular, for a;'' = t/)'^, and for the appropri- 
ate function U , one recovers the Fierz and Pauli 
(1939) mass term. 

De Rham et al (2010) show that a two parame- 
ter family of massive gravities has some good fea- 
tures (ghost-free in the so-called decoupling limit). 
The family is described by the potential (Nieuwen- 
huizen 2011, Hassan, Rosen 2011, Koyama et al 
2011) 

[/ = ^ \a^b + C2 ^ AyiAsAc + C3A0A1A2A3, 

. 

where the sums are over all all-distinct pairs and 
triples of indices, and A^i are the four eigenvalues 
of the matrix 




We call this two-parameter family of theories 
FP2. 
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B. Static spherically symmetric field 

Starting from the unitary gauge x'^ — cj)'^, the 
asymptotically flat metric can be written as (Ap- 
pendix E) 

ds^ = e^df ^ e^dR^ ~ R^et'dn^. (Bl) 

We then change the radial coordinate to transform 
the metric into the usual form 

ds'^ = e^df - e^dr^ - r^dfl'^. (B2) 

The scalar fields in the new coordinates become 



(B3) 



where is the unit radial vector. 

As independent equations we use two Einstein 
equations and the stress-energy conservation 

G"„ = r-2(i _ e-^) + r-^e-^X' = T^, (B4) 

Gl^r-^{l-e-^)-r-^e-^iy' ^T^ (B5) 

7^1 ' = I'^'iT^ - Tl) + -(T| - Tl)^ (B6) 
2 r 

where prime denotes the r-derivative. 

The stress-energy tensor of FP2 is derived from 

ea. (|M)) : 



T° = -m2(Ai + 2A2 + c2(2AiA2 + A^) + cgAiA^), 

(B7) 

Tl = -m2(Ao + 2A2 + c2(2AoA2 + A^) + cgAoA^), 

(B8) 

= -rr?{\o + Ai + A2 + (B9) 
C2(AoAi -|- A0A2 -l- A1A2) -f C3A0A1A2), 



with 
and 



C2 = 1 + C2, C3 = C2 -I- C3, (BIO) 



Ao = l-e-''/^ (Bll) 

Ai = 1 - ^ ^ _ e-^'^{l - A2 - rA^), (B12) 

A2 = l-e-^/^ (B13) 

We should remark that e"^/^ in (|B3IB13p and 
ij-^ (|B12|1 are defined as square roots of 
and and can therefore be negative, because 
the square root function must be allowed to switch 
to negative branch at zero as continuity dictates. 



We will ultimately work in terms of A2 instead of 
and w = e"'*'/^, which automatically selects the 
right sign of the square root. 

Using expressions (|B7IB8IB9[) and expression 
(|B11IB12P in eqs. (|B4IB5IB6p gives the following 
system of three equations for three unknowns 
A, and A2: 



1 - e-"- + re^-'A' = -TO^r''(2A2 + C2A^ + 
(1 + 2C2A2 -I- C3A2)(1 - e-V2(i _ A2 - rA'2))) 



(B14) 



1 - e"^ - re-^v' = -m^T^{2\2 + C2A^ + 

(l + 2c2A2+C3A2)(l-e-''/2)) (B15) 

rv'(\ + 2C2A2 + CaA^) ^ 4(e^/2 _ ^^(^ ^ ^^y^ ^ 

(c2+C3A2)(l-e-'^/2)). (B16) 



After linearizing in A, /i, equations (jB14| 
IB15| IB16P can be solved exactly: 



-ce 



^(-+m), (B17) 



2m2 



r 

„ , 1 ra rr? , 



manifesting the vDVZ discontinuity. 

A Vainshtein type argument can be used at 
this point. Neglect the graviton mass in eqs. 
(|B14IB15P to recover the Einstein theory. Then 
substitute the linearized Schwarzchild solution 
V = —A = —Tg/r into the Bianchi identity (jB16l) . 
also linearized in v and A, but exact in A2, and 
solve for A2. One finds 



C3A2 



1. 



(B18) 



which implies C3 > 0. 

Without hnearizing, equations (|Bl4l [BTsI IBIHI) 
cannot be solved as written. Although we do have 
three equations for the three unknowns, it is seen 
that the only derivatives of the unknowns in these 
equations are A' and Aj in eq. (|B14[) . v' in eq. (jBlSp , 
and v' in ea. (|B16p . 

Equating the two expressions for u\ we derive 
an algebraic relation between r, A, A2. Then one 
can select any two of the three unknowns v, A, A2 
and derive a system of two first-order differential 
equations for the two selected unknowns. 
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In practical, we have used 



(B19) 



as the two unknowns. Then the algebraic relation 
can be written as 



1 - e 



-u/2 



F 
G 



(B20) 



where 



F = (1 + 2c22/ + C32/')(l -w^+ m\^{2y + cay')) 

-A{l + C2y)w{l-w), (B21) 

G = ^C2+C3y)w{l-w)-m'^r^{l + 2c2y + C3y^)^. 

(B22) 

and the final form of the two equations for the two 
unknowns follows from (|B141IB16( ) and (|B20p : 

aiw' + biy' = di (B23) 

a2w' + b2y' = d2 (B24) 

Here (|B23l) is just (|BT41) in the new notation, 
so that 

ai = 2rw (B25) 
bi = -m'^r^{l + 2c2y + C3y^)w (B26) 

di = l-w'^+ ■m'^r^{2y + C2y^ + 

{l + 2c2y + c-iy^){l-w + wy)) (B27) 

To calculate the coefficients 02, 62, (i2, we first 
write eq. (jB16p as 



where 



azy'=/3(l-e~''/^) + 7 (B28) 

a = rw[l + 2c2V + c^y"^) (B29) 

/3 = 4(l-w)(c2+C32/) (B30) 

7 = 4(l-ii;)(l + C2y) (B31) 



We then substitute 1/ and i^' , as given by (|B20p , 
into (jB28p . We will need the derivat ives: 

F' ^ Fiw' + F2y' + F3 (B32) 

G' = Giu;' + G22;' + G3 (B33) 



^ The results of the following calculation are to be used by a 
computer. Therefore, we organize the calculation so as to 
minimize the chance of human error, rather than to derive 
formulas that look nice to the human eye. 



where 

Fi = -2u;(l + 2c2y + c^^y^) - 4(1 + C2y)(l - 2w) 

(B34) 

F2 = 2(c2 + C3y)(l -w'^ + m^r^{2y + C2y^)) + 
2(1 + 2c2y + C3t/2)mV2(l + C2y) - 4c2(u; - w^) 

(B35) 

F3 = 2m\{l + 2c2y + C3y'')[2y + cay^) (B36) 

Gi =4(c2+C3t/)(l-2u;) (B37) 

G2 = 4c3(w-w^)-4mV2(l+2c2y+C3y2)(c2+c3?;) 

(B38) 

G3 = -2mV(l + 2c2y + c^y^f (B39) 
Then (jB28l) and (|B201) give 

a2 = 2a(GFi - FGi) (B40) 

62 = 2a(GF2 - FG2) (B41) 
d2 = -2a(GF3-FG3) + (G-F)(/3F+7G) (B42) 

For black holes, the r-derivatives in equations 
(|B23|B24p are converted to the p-derivatives by 



wd/dr — d/dp . 



(B43) 



C. Static spherically symmetric field with 
matter 

To generalize the formalism of Appendix B in 
the presence of a perfect fluid one adds 



= diag(e, -p, -p, -p) 



(CI) 



to the r.h.s. of Einstein equations which changes 
(|B4|B5P and (|B14|B15P accordingly. Besides, 
there are now two copies of stress-energy con- 
servation (|B6p . one is for the mass term (|B16[) 
and the other one 



(C2) 



describing the hydrostatic equilibrium. The equa- 
tion of state p(e) and (jC2[) complete the system 
(|B14I IB15I IB16p . to be solved for five unknowns 
V, A, A2, e and p. 

We still have two equations with v' as the only 
derivative. This allows to express v and v' in terms 
of A, A2, and p: 



l-e 



-u/2 



F/G 



(C3) 
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with 

F ^ F + (1 + 2c2y + C3y^)pr2 , (C4) 

and 

z/' = /3F/(aG) +7/a, (C5) 

where F, G, a, /? and 7 are given in Appendix B. 

The rest of the calculation proceeds with minor 
modifications to yield 

aiw' + biy' = di , (C6) 
a2w' + h2y' + i2v' = ^^2 , (C7) 

which can be integrated together with (|C2[) . Here 

the new coefficients are given by ai = ai , 61 = 61 , 
a2 = 02 and 

Ji = (ii - er2 , (C8) 

02 = 2a(GFi -J^Gi), (C9) 

62 = 2a(GF2 - FG2) , (Cf 0) 

d2 = -2a{GF:i~FGz) (Cff) 
+ (G-F)(/3F + 7G), 

e2 = 2ar2(l + 2c2y + C3y2)G, (Cf2) 

with 

i^i = , (Cf 3) 

F2 = F2 + 2(c2 + C32/)pr2 , (Cf4) 

F3^Fs + 2{l + 2c2y + C3y^)pr. (Cf5) 

D. The near- horizon solution 

We want to find analytic expressions for the 
near-horizon solution in the small-graviton-mass 
limit. We first find the eigenvalues Xa- Since Ai 
and A2 remain finite at the horizon, they can be 
replaced by their zero-graviton-mass limit. Substi- 
tuting = = f — Tg/r in the Bianchi identity 
(|Bf 6P and solving for A2(r) gives, at r = rjf|: 

A2 - -2 - C2/C3 + ^4 + 4/cl - f /C3 , (Df ) 
rg^'2 = -A2 - 2c2/y/4ci + ci-C3, (D2) 

here and below the prime denotes the p-derivative, 
with dr/dp = w — e^^^^. The same substitution 
= f — Vg/r in (|B12[) gives, 'dt r — Vg, 

Ai = f + rgX'^ . (D3) 

*The larger root was chosen for A2 because in the "Vainshtein 
region" A2 = l/^/c3 which is larger than both roots. 



Now we can find the desired solution. Due to 
finiteness of Ai and A2, the contribution of the 
mass term to {00} component of Einstein equa- 
tions (jBf 4p can be neglected. Thus we have 

2u/ = (f - w2)/r ~ f/r, (D4) 

which gives w ~ p/2r^in and 

r~rmi„-t-pV(4rmi„). (D5) 

Knowing w and y — X2, the redshift z = e^"/^— 
f = — Aq can be calculated from (|B20[) and is of 
the form 

z ~ C4/{p/rnun + C^m^^) , (D6) 

where G4 and G5 are positive dimensionless num- 
bers for the parameter range of interest. 

Finally, to derive an expression for the curva- 
ture, we note that the near-horizon T^f is domi- 
nated by the Aq = —z contributions: 

Ti~m2z(f + 2c2A2+C3Ai), (D7) 
T| ~ m^z{l + C2(Ai + A2) + C3A1A2) . (D8) 

Requiring that the leading order given by 
(|D7p vanishes at the A2 value given by (|Dip , gives 

C3 = cl (D9) 

And for C3 = C2, 

r| ~TO2^(f -HC2Al)(f +C2A2) (DfO) 

also vanishes at the A2 value given by (jPf p . 

E. The asymptotically flat metric 

To define the graviton's potential in Lorentz in- 
variant theories of massive gravity one has to in- 
troduce a flat reference metric rj^^, . In the unitary 
gauge (x^ = (/)^) in which we work, the appearance 
of ry^ii/ breaks the reparametrization invariance of 
general relativity. Nevertheless the theory is still 
invariant under simultaneous transformations of 
the physical metric and the reference metric, and 
therefore, it is convenient (and always possible) to 
take the latter to be of the form: 

rj^^dx>'dx'' = dt^ ~ dr^ - r^dV? , (El) 
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when dealing with a sphericahy symmetric setup. 
The most general static metric which respects the 
symmetries of the problem is 



ds^ = e^df - e^dr^ - 2£,dtdr - 



(E2) 



where i/, cr, ^, and fj, are arbitrary functions of r. 

Here we show that the requirement of asymp- 
totic flatness, in the sense that at r — )■ oo, g^^, ap- 
proaches the flat metric as deflned in (jEip . forces 
^ to be identically zero. 

The key observation, first made by Salam and 
Strathdee (1977), is that the identity gtrRtt — 
guRtr = which holds for the Ricci tensor of (IE2I) . 
imposes a purely algebraic constraint on the com- 
ponents of the metric 



gtrTtt — guTtr — , 



(E3) 



with being the effective stress-energy tensor 
of the mass term. 



For the metric (jE2p . and arbitrary potential 
that is a function of the eigenvalues of the matrix 
in (jA4p one has Ttr = ngtn with the coeflicient 
K which is non-singular at gtr = 0. It follows that 
eq. (|E3|) divides the solutions into two categories 

(i) g^^ = 0, 

(n) Ttt = K.gtt- 

Now consider large values of r, where the met- 
ric perturbations h^^ = giiv — Vp.i' is small by the 
asymptotic flatness assumption. Then, for any 
non-linear completion of the Fierz-Pauli mass term 
/i^, we have 



Ttr^-^m^htr{l+Oih)), (E4) 

giving K — — at large r. This excludes the 
branch (ii) as it leads to a finite value for Tu- 

The branch (ii) is the Schwarzschild-de Sitter 
solution of Salam, Strathdee (1977) and Koyama 
et al (2011). 
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